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ABSTRACT 

The matrix sine- Gordon theory, a matrix generalization of the well-known sine- Gordon 
theory, is studied. In particular, the ^-generalization where fields take value in SU(2) de- 
scribes integrable deformations of conformal field theory corresponding to the coset SU(2) x 
SU{2)/SU{2). Various classical aspects of the matrix sine-Gordon theory are addressed. 
We find exact solutions, solitons and breathers which generalize those of the sine-Gordon 
theory with internal degrees of freedom, by applying the Zakharov-Shabat dressing method 
and explain their physical properties. Infinite current conservation laws and the Backlund 
transformation of the theory are obtained from the zero curvature formalism of the equation 
of motion. From the Backlund transformation, we also derive exact solutions as well as a 
nonlinear superposition principle by making use of the Bianchi's permutability theorem. 
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1 Introduction 



The sine-Gordon theory is the most well-known example of relativistic integrable field the- 
ories in 1+1 dimensions. Exact solutions, topological solitons and breathers, are known 
and various applications of these solutions have been made in the study of a wide range of 
physical systems. However, in many cases, such an application has been made after a trun- 
cation of the physical systems, i.e. suppressing all degrees of freedom except one scalar field, 
so that the sine-Gordon theory becomes an effective description of the reduced system. In 
this regard, it is desirable to have generalizations of the sine-Gordon theory with additional 
degrees of freedom whereas the integrability of the theory is maintained to provide exact 
solutions. One such example is the so-called complex sine-Gordon theory which appeared 
in the description of the relativistic vortex motion in a superfluid[Q, in a treatment of 0(4) 
nonlinear sigma model @ and more recently in the context of conformal field theory 0. In 
the complex sine-Gordon theory, the scalar field is complex valued so that the phase factor 
becomes an additional physical degree. Other types of nonabelian generalization of the sine- 
Gordon theory, carrying extra internal degrees, have been considered recently in the context 
of integrable deformation of the coset conformal field theory [Q, particularly in the case of 
the critical Ising model deformed by $(2,1) and $(3,1) operators 0. In fact, to each classical 
Lie algebras, there exist nonabelian generalizations of the sine-Gordon theory which admit a 
positive definite kinetic energy when certain criteria of the sl(2) embedding are metH. The 
complex sine-Gordon theory and the critical Ising model deformed by $(2,1) operator arise 
as first two examples in these generalizations which correspond to the Lie algebra B2 and 
A 3 respectively. 

In this Paper, we study in detail the ^-generalization of the sine-Gordon theory, which 
we call "the matrix sine-Gordon theory", at the classical level.0 At the quantum level, this 
theory may be regarded as an integrable deformation of a minimal model M corresponding 
to the coset SUn(2) x SUn(2)/ SU2n{2), in particular the deformation of the critical Ising 
model by the energy operator $(2,1) if the level N = 1. We demonstrate the classical inte- 
grability of the theory by deriving infinite current conservation laws from the zero curvature 
formalism of the equation of motion. The Backlund transformation is also obtained from the 
zero curvature formalism, and using the Bianchi's theorem of permutability|J a nonlinear 
superposition principle is derived for the solutions of the matrix sine-Gordon theory. Exact 
solutions, solitons and breathers which generalize those of the sine-Gordon theory, are ob- 
tained by making use of the Zakharov-Shabat dressing method|l(|. These solutions carry 
internal degrees of freedom which affect the scattering process. Two soliton solutions for the 
soliton(antisoliton) - soliton(antisoliton) scattering is given explicitly in terms of parameters 
u and 9 which describe the relative velocity of solitons and their relative internal directions 
respectively. It is shown that the nonabelian effect which is controlled by the parameter 9 
makes solitons less repulsive for larger 9 while the scattering time of far removed solitons 

3 The name "matrix sine-Gordon" has been also used for a different model on the group 
SO{n) in §. 
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does not depend on 9. An explicit form of nonabelian breather solution is also given in 
terms of both energy configuration and internal components with parameters v and 9. The 
breather solution is a bound state of soliton and antisoliton which oscillates in time with 
angular frequency 2^—kv/ yl + v 2 . As 9 increases, the breathing mode of the potential 
energy configuration diminishes and at 9 = ir, it becomes completely breathless! However, 
the internal components still oscillate at 9 = it which keeps the static potential energy con- 
figuration. Alternative derivation of one and two soliton solutions are also given using the 
Backlund transformation and the nonlinear superposition principle. 

The plan of the Paper is the following; in Sec. 2, we define the model in terms of the 
gauged Wess-Zumino-Witten functional and address issues on the gauge invariance and the 
gauge fixing. We present the zero curvature formalism of the equation of motion for the 
model and from which, derive infinite current conservation laws thereby demonstrating the 
integrability of the model. In Sec. 3, exact solutions, N solitons and breathers, are derived 
rigorously by making use of the dressing method. A pictorial description of these solutions is 
given and their physical meaning is discussed. Sec. 4 deals with the Backlund transformation 
and the nonlinear superposition principle and Sec. 5 is a discussion. 



2 The Model 

We introduce the matrix sine-Gordon theory from the context of conformal field theory and 
its integrable deformations. The action principle for the G/H-coset conformal field theory 
may be given in terms of the gauged Wess-Zumino-Witten(WZW) functional [11], which in 
light-cone variables is 

S(jg, A, A) = S wzw {g) + ±- [ T^-Adgg- 1 + Ag^dg + AgAg' 1 - AA) (1) 

Z7T J 

and Swzw(g) is the action of group G WZW model 

Swzw(g) = ~ [ Tr g- x dgg- x dg - -L [ Tr g^dg A g~ x dg A g^dg . (2) 

47T JS 127T JB 

g is an extension of a map g : X — > G to a three-dimensional manifold B with boundary 
g\dB = g, and the connection fields A, A gauge the anomaly free subgroup H of G. Here, 
we take the diagonal embedding of H in Gl x Gr, where Gl and Gr denote left and right 
group actions by multiplication (g — > giggR 1 ), so that Eq.(l) is invariant under the vector 
gauge transformation; g — > hgh^ 1 with h : £ — »H. In particular, the minimal unitary series in 
conformal field theory arise from the restriction to the coset, (SU(2) L x SU (2) N ) / SU (2) L+N , 
where integers L, N denote the level of the Kac- Moody algebra|8]]. In this case, the full theory 
is given formally by functional integrals, 

J [d gi ] [dg 2 ] [dA] [dA] exp il ( gi , g 2 , A, A) (3) 
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where 

I { gi , g 2 , A, A) = LS wzw ( gi , A, A) + NS wzw (g 2 , A, A) (4) 

and A, A gauge simultaneously the diagonal subgroups of SU(2) x SU{2). 

The matrix sine-Gordon theory is defined as a massive deformation of the minimal series 
with L = N which, at the action level, is adding a potential term to the action, 

Nk f 

I(9i,92, A, A, re) = I Q (#i, 02, A, A)- — / Tr(g± 1 g 2 + g 2 1 9i), (5) 

Zix J 

where re is a coupling constant .0 Note that the potential term is invariant under the similarity 
transform; g\ — ► sgis~ l , g 2 — > sg 2 s~ 1 so that the vector gauge invariance of the action is 
maintained. In the convention of coset conformal field theory, the potential term transforms 
at the classical level as (doublet, singlet) so that it corresponds to the integrable perturbation 
of minimal series by the operator ^(2,1) 123 ■ ^ n ^ ms P a P er 5 we focus only on the classical 
aspect of the theory so that the level N is irrelevant. In order to understand the vacuum 
structure, we parameterize 1 g 2 by 

3 

9i 1 92 = exp(icj)&) , a = Y^ a i a i ( 6 ) 

i=i 

where Oi are Pauli matrices and the coefficients are normalized to one, X)f=i a i a i = 1- 
Then, the potential becomes 

xr Nk x _ x 2Nk 

V = —Tr{g 1 g2 + 9 2 9i) = cos0 . (7) 

Z7T 7T 

If the coupling constant re < 0, V possesses degenerate vacua at cf) = 2mr for integer n so 
that g\ X g<i = 1 or for any arbitrary g\ = g 2 valued in SU(2). Note that a specific vacuum 
is characterized only by the integer n independently of eij. The degeneracy of the vacuum 
allows soliton solutions which interpolate different vacuua. The explicit solutions will be 
found in Sec. 3 and Sec. 4. The topological soliton numbers are defined by the difference 
An = rii — n 2 of integer values of two interpolating vacuua. If re > 0, degenerate vacuua 
occur at <\> = (2n + \)ix for integer n so that g\ X g 2 = —1, or for any arbitrary g\ = —g 2 . 
From now on, we will restrict ourselves to the re < case only, re > case will be discussed 
in Sec. 5. 

The classical equation of motion arising from the action Eq. (5) is 

[ d + g x 1 3g 1 + gl x Ag x , 3 + A}- n(g x 1 g 2 - g 2 1 g 1 ) = 

[ d + g 2 1 dg 2 + g 2 1 Ag 2 , 8 + A } + K{g x l g 2 - g 2 l 9l ) = (8) 

whereas variations of the action with respect to A and A give rise to the constraint equation, 

-dgigx 1 + g l Ag 1 - 1 -dg 2 g 2 - 1 + g 2 Ag 2 1 -2A = 

g^dgi + 9i 1 Ag 1 + g 2 l dg 2 + g 2 l Ag 2 - 2A = . (9) 

4 Here, we asuume gi and (72 to take values in SU(2). The U(2) case has been considered 
in § . 
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The constraint equation in Eq.(9), when combined with the equation of motion in Eq.(8), 
results in the flatness condition of A and A, 

dA-BA+[A,A] = 0, (10) 

which reflects the vector gauge invariance of the action. In the following, we consider two 
types of different gauge fixing. Assume that the underlying manifold £ is the flat two- 
dimensional Minkowski space R 1+ . Then the flatness of A and A allows us to choose a 
"nonlocal gauge"; f\ 

A = A = 0. (11) 
The equation of motion in the nonlocal gauge becomes 

%i~ 1<9 #i) + K{gi 1 g 2 - g 2 1 9i) = 

B(g^dg 2 ) - K(g^g 2 - g 2 1 g 1 ) = (12) 

whereas the constraint equation becomes 

dgigi 1 + dg 2 g 2 x = , g^ 1 dg l + g 2 1 dg 2 = . (13) 

In the abelian limit, where g\ = expi<fiiO~i, g 2 = expi<j) 2 ai, the constraint equation may 
be solved locally by <pi = <ft — — 02 so that the equation of motion in terms of <fi becomes 
precisely the sine-Gordon equation. However, for g±,g 2 valued in SU(2), the constraint 
equation in general can not be solved locally. Consequently, in the nonlocal gauge, a local 
parametrization solving the constraint is not possible as in the abelian case.[] 



Nevertheless, there exists another type of gauge fixing, so-called "the unitary gauge" [[13 
which allows a local parametrization solving the constraint in the following sense; for any 
given gi and g 2 , one may bring g 2 into a form g 2 = exp(i0o"3) via the similarity transform; 
9i —> s 9i s l ,92 sg^ 1 for some s. Thus, the scalar function cf) parameterizes the equiv- 
alence classes of g 2 with the equivalence relation given by the similarity transform. The 
remaining U(l) gauge symmetry which leaves g 2 invariant may be used to fix g\ to give the 
unitary gauge; 



_ ( ue~ l<t> Wl~uu*e^\ _ fe^ \ . . 

In this gauge, the constraint equation can be solved explicitly for A and A in terms of u and 



5 For general E, such a nonlocal gauge is not always possible because holonomies of 
A, A could be nontrivial. 

6 Even though local parametrization is not possible, exact solutions can be constructed 



explicitly as in Sec. 3 and Sec. 4. 
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where 



b = 



-21 



^du idu + idu* 



2\J\ — uu* 8vl — mm* sin 

(e 2i(t> u* - u)du* - (e~ 2i(t, u - m*)<9m (m - u*) (du + 9m*; 



and 



4(1 -mm*) 16(1 — mm*) sin 

du* ie 2l(t> du* — ie~ 2t<i> du 



(16) 



4\/l — mm* sin 2 8 vT^ml* sin 2 

( e ii<t> u * _ _ {e~ 2i< t>u - u*)du (u - u*)(du + du*) 

a = —lOd) ; r 1 ; r 5——. I 17) 

Y 4(1 -mm*) 16(1 — mm*) sin y ' 

Then the flatness condition Eq.(lO) resolves into the equations of motion for u, u* and 0. 
These equations are related to those of the nonlocal gauge by the following association; if we 
solve Eq.(lO) in terms of holonomies, A = h~ l dh and A = h~ 1 dh, then g 1? g 2 of the nonlocal 
gauge are related to those of the unitary gauge by 

g? = hg^h- 1 , g? = hg u 2 h~ l . (18) 

For the rest of the Paper, we will restrict ourselves only to the nonlocal gauge. Translation of 
subsequent results into the unitary gauge can be readily made by the association in Eq.(18). 

In order to understand the integrability of the matrix sine-Gordon theory, we consider 
the linear 4x4 matrix equations with a spectral parameter A, 

Li(A)tf = (d + U — AT)* = 

L 2 (A)^ = (d + A+W 1 )* = (19) 

A 

U = Q~ l dQ + Q~ X AQ , V x = Q~ X TQ , T = ikE , f = iS (20) 



where 
and 



with each entries being 2x2 matrices. The matrix sine-Gordon equation arises precisely 
as an integrability condition, [Xi(A) , -^(A)] = 0, of the linear equation for any A. The 
advantage of the linear equation with a spectral parameter A is that it allows a systematic 
way to construct infinite conserved currents. In addition, exact solutions can be obtained 
from the linear equation which we consider in Sec. 3. In order to find conserved currents, we 
solve the linear equation iteratively by setting 

oo 

$ = ^ exp(-AT^) = X ~ m $m ; $0 = 1 (22) 

m=0 
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so that the m-th order equation in the nonlocal gauge is 

<9$ m + Uq<& — [T, <& m+ i] = 
<9$ m + Vi$ m _i = 
rri _ 1 ( Q^dgi - g 2 1 dg 2 \ , , 

u ° ~ 2 1 o - ^r 1 ^! J • (2d) 

In accordance with the initial value $ — 1, can be appropriately parametrized by 

p 2 m \ _ _ / 



•~ s { o s j • u: +1 rj""^ - < 24 > 

This brings Eq.(23) into a component form, 

<9P m+ i + igi 1 g 2 S m = 
dSm+i + ig^giPm = 
<9P m + gi 1 dg 1 P m = in(P m+1 - S m+1 ) 

dS m + g 2 1 dg 2 S m = iK(S m+1 - P m +i) , (25) 
which we solve iteratively with initial values P = S = 1, 

P m+1 - S m+ i = ^—\d(P m - Sm) + g^dg^Pm + S m )\ 

P m +i + Sm+i = -J dz(ig^ 1 g 2 S m + ig 2 1 g 1 P m ) - J dzg^ 1 dg 1 (P m+1 - S m+1 ) . (26) 
In particular, 

Pi = ^~9i l dgi - ^ / dz^dgif - % -J dz(g^g 2 + ^ V) 

Si = -^ 1 dg l -^-Jdz(g^dg 1 ) 2 - i -Jdz(g^g 2 + g 2 1 g 1 ). (27) 

With iterative solutions of the linear equation, we find that the consistency of Eq.(25), 
dd P m = BdP m and ddS m = BdS m , gives rise to two sets of current conservation laws; 

<9J« + 9M, = 



'm+2 
r(2) 



where 



5j£> + 9Ji 2 | 2 = 0; m>0 (28) 
•^m+2 = 9P m+ i = —g 1 1 dg 1 P m+ i + in{S m+ \ — P m+ \) 

Jm = w^giPm 

Jm+2 = 9S m+ i = g 1 1 dgiS m+ i + in(P m+ i — S m+ \) . (29) 
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The subscript m of the current J m denotes the conformal spin in the massless limit, or it 
simply counts the order of the derivatives. In particular, the m = case is the energy- 
momentum conservation, 

BT ± + d<S)± = 0, (30) 

where 

T + = «e(jf + J 2 (2) ) = (g^d gi ) 2 

+ = i K (jU + jW) = - K ( gI lg 2 + g^ gi ) (31) 

while the other half of the conserved currents are 

T _ = iK (jV - jf) = d(g^d 9l ) 

Q_ = iK{jW-J®) = -K( a ?g 2 -g?g x ). (32) 

It is interesting to observe that T__ in the abelian limit becomes T_ = d 2 <ft which is precisely 
the term added to improve the energy- momentum tensor in the Feigin-Fuchs construction Jyj] . 
Another type of conserved currents arises from the invariance of the matrix sine-Gordon 
theory under the parity transform, 

z <-> z and g x <-> g^ 1 , g 2 g^ 1 ■ (33) 

This leads to the parity conjugate pair of conserved currents which, together with currents 
in Eq.(29), constitute a complete set of conserved currents of the matrix sine-Gordon theory. 
For example, the parity conjugate of the energy-momentum is 

f + = {dg l g^)\Q + = -^g l g 2 l +g 2 g^) (34) 

and 

T_ = -d{dg x g x l ) , 0_ = -K{ 9l g 2 l - g 2 g^) . (35) 



3 Dressing Method and Soliton Solutions 

In this section, we give a detailed account of the derivation of soliton solutions. We follow the 
dressing method of Zakharov and Shabat|10| and obtain nontrivial soliton solutions from the 



trivial one by employing the Riemann problem technique with zeros ||15||. In Sec. 4, we give 
an alternative method based on the Backlund transformation and obtain soliton solutions 
by direct integration. We first give a brief review on the dressing method. For later purpose, 
we rewrite the linear equation in the nonlocal gauge by making a similarity transform of 
Eq.(19) by the matrix Q, 

« = «"' S ^G -0 (36) 



8 



such that 
where 



(d + U' - AT>' = , (8 + = (37) 



tf' = QVQ- 1 
U' = QUoQ' 1 

and 

T , = QTQ- 1 = <«(J _°J , f' = QfQ- 1 = t(J (39) 

In the following, we drop the prime for convenience without causing any confusion. The 
dressing method is a systematic way to obtain nontrivial solutions from a trivial one. In our 
case, we take the vacuum as a trivial solution of Eq.(37), 

gi = g 2 = l and ^ = ^° = exp(AT^ - A _1 T z) . (40) 

Let T be a closed contour or a contour extending to infinity on the complex plane of the 
parameter A. Consider the matrix function $+(-2, z, A) which is analytic with n simple poles 
/i!,...,/i n inside T and $_(^,z, A) analytic with n simple zeros Ai,...,A n outside T. We 
assume that none of these zeros lies on the contour T and \[/ (\[/ )~ 1 = ($_) _1 $ + = 1 for 
A 7^ fii, \i ; i = l,...,n. We normalize $+, $_ by $ + |a=oo = $-^=00 = 1- Differentiating 
\E f0 \l/ ~ 1 = ($_) _1 $ + = 1 with respect to z and z, one can easily see that 

- i$+T &Z 1 = d^^Z 1 - j$-T $l\ (41) 
A A 

Since $+($_) is analytic inside (outside) T, we find that the matrix functions Uo and V\, 
defined by 

U = -d®®" 1 - $AT$ _1 + AT ; VI = -A0$$ _1 + ^T^ 1 (42) 

where $ = $ + or <3>_ depending on the region, become independent of A. Also, ^ = 
satisfies the linear equation; 

(d + U - AT) 1 !' = , (9 + T Vi)* = . (43) 

A 

The identification \& = \I> and Uo, V\ with respect to Uo, V\ in Eq.(38), C/o and V\ then 
provide nontrivial n-soliton solutions. 



g^dgi 
g^dgi 



l(9i 92 

2 UrV 



92 1 9i 
92 X 9\ 



92 1 9i " 
-91*92 



9i l 92 



(38) 
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In making such an identification, the specific form of Uq and V\ imposes restrictions on 
\l/ which in certain cases may be solved algebraically. For example, the anti-unitarity of U 
and V\ imposes restrictions on \1> and <£> which may be complied with 

tft( A ) = ^-i(A*) , $t(A) = $- x (A*) . (44) 

These are not the most general expression giving the anti-unitary Uq and Vi, however we 
assume Eq.(44) since they suffice for our purpose of deriving soliton solutions. In order to 
construct the matrix function <3> for the soliton solutions, we take the ansatze for $ and 

M Aa M -na 

* = i + ECt^V) > ^ = 1 + (45) 

a=l A ~ Ua a=l A ~ ^ 

where the matrix functions A a (z, z), B a (z, z) are to be determined. Since the identity; 
<3><3> _1 = 1 and Eq.(42) should hold for any A, they require respectively algebraic and differ- 
ential relations among A a and B a . These relations can be obtained through the evaluation of 
residues of both equations at A = ji m v a . For instance, the residues of the equation = 1 

gives rise to 

M d/3 M Af3 

A a + A a J2( ) = , B a + J2( )B a = 0, (46) 

{3=1 v a - p =l Ha ~ V/3 

while those of Eq.(42) lead to 

M r>P M a/3 

A a D h2 (u a )[l + £( )] = , [1 + E( )]DiMB a = (47) 

/3=1 u a - HP /3 =1 Ha - Vp 

where 

Di{\) = d - XT , D 2 {\) = B + \f. (48) 

A 

In order to solve Eqs.(46) and (47), we assume that Afj = s l a t J a , Bf^ = n^m^ where 
s ai^ai n ai m a are ^ wo ^ w0 matrices with i — 1,2. Then, Eqs.(46) and (47) changes into 

2 M n 3 2 M J +i 

E W + E -^-] = o , E^' 1 + E ^-H = o (49) 

j=l 0=1 U <* ~ VP i= l p=i Hs ~ Vp 

and 

Di, 2 (Ha)np = , tpD h2 (u a ) = , (50) 

where we understand tpd and tpd as — dtp and —dtp. Note that n a and t a can be solved in 
terms of arbitrary constant vectors n a and i a , 

< = Et*°wrx , t i a = iimvaVrti , (si) 

3=1 3=1 
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while m a and s a can be obtained in terms of t a and n a by solving the linear algebraic 
equation (49) such that 

M M 

< = -E > 4 = E > ( 52 ) 
0=1 0=1 

where 

"^E-^-. (53) 

V -1 is defined by Y!p=i VapVp-y — <W1 where 1 is the 2x2 unit matrix. The unitarity 
condition Eq.(44) requires that 

Ha = S a and A ajf — B a ( or # = <, 4 = < )■ (54) 

Consequently, 

€ = < , (Va) f = - V • (55) 

Further specification of t a and n Q arises from the identification; U = U , V\ = Vi, in Eq.(42). 
Since Eq.(42) holds for any A, we combine Eq.(42) and (45) and take the A — > 00 limit to 
obtain ^ : ^ ^ 

TT — A ar T _l_ T n a \ — ( S ata n a m a ~~ S a^a n a m a\ ( KR\ 

~ U - 2^\ A 1 + 1 a ) - lK l^ [ q 2 f l _ 2 1 _„2 f 2 _ 2 2 J • l^OJ 

Note that the (block)-diagonal part vanishes identically due to the equality 

+ « = "J^" 1 W« - rtjy- 1 )^ = (57) 
which agrees with Uq. The off-diagonal part gives rise to 

g?d 9l = -i«(-4*a + «) = 2«E#« = 2i«E^> (58) 

or 

sr 1 ^ = -<«(4*i - = -»« E = -2™ E ^ , (59) 

where Z Qi g and Z a( g are defined by 



J aP 



Vp — Ha Hp — He 



jyg = nj(nj)- 1 = exp(2A /3 )nJ(n2)- 1 = exp(2A /3 )iV /3 

%z 

A/3 = infipz (60) 

Hp 
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and we have used the unitarity condition Eq.(54). The last step in Eq.(58) (similarly Eq.(59)) 

1(3 V a f3 tfi 



can be checked easily by using K\ = J2b V a 3"^8 so that 



-n. 



(61) 



Thus, the identification with Uq through Eqs.(58) and (59) imposes restrictions on Z a p such 
that 

E^-E^E^ 1 - (62) 

a,0 a,3 a, 8 

On the other hand, the V\ part in the A — > limit of the linear equation, gives rise to 

V 1 = $f $-^=0 or $| A=0 f = V^\ x=0 . 



In components, they are 



9i 92 = -($12 + $22) ($12 - $22) 



-1 



(63) 



(64) 



or 



9l 1 92 = -($11 + $2l)($21 - $ll) 1 



(65) 



where <3>jj denote 2x2 block components of $|a=o- We will see below that the two expressions 
for 9\ X 92 in Eqs.(64) and (65) are indeed equivalent when the condition Eq.(62) is satisfied. 
If we define 

_ l+NiN p 

*a8 = 



1 - aWa£ ' 



Eqs.(64) and (65) become 



9^92 = -($12 + $22) ($12 - $22) 1 



i-E 



i-E 



i-E( 1 + ^)( y_1 W i-E( 1 -^)(^" 1 W 

q,/3 a,/3 



and 



0i £2 



($11 + $2l)($21 - $ll) 1 



i + E 



-4 + 4)£ 



1 - EC 1 + + E(! - JV/jX^W^ 



(66) 



(67) 



(68) 



These are M-soliton solutions of the matrix sine-Gordon theory. In the following, we give an 
explicit expression for M=l and 2. 
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M=l; 1-soliton 



For M=l, we have 

Z=-tf = (69) 

and Eq.(62) for M=l, Z~ x = -Z~\ can be solved either by N = -N' 1 or N = -Nl The 
former case results in only a trivial solution while the latter case, N = exp(2A)iV = — = 
— exp(2A*)N\ requires A real and N anti-hermitian. Thus, \i = i8 is pure imaginary and 
A = —k,5z — z/5. We parameterize the anti-hermitian matrix N by N = i exp(^)aj(jj where 
(Tj are Pauli matrices and the repeated index i denotes summation from i=l to 3. 77, Oj are 
arbitrary real constants with a normalization = 1. Then, from Eqs.(67) and (58) we 
obtain the 1-soliton solution given by 

= (r^) 2 - (- ta ° h ( 2A+ ") +io " r ' ^h(^) ) ■ (70) 

and 

^r 1 ^! = 2mZ~ 1 = 2m<5 — — -a^ . (71) 

cosh(2A + rj) 

Combining Eqs.(70) and (71), we could solve for gi and g<i- They agree with the explicit form 
given in Eq.(131) of Sec. 4 which is derived directly from the Backlund transform. Physical 
meaning of parameters in the soliton solution is the following; parameter rj depends on the 
choice of origin of space and time. We choose the origin to set them to zero and introduce 
the space and time coordinate by t = z + z, x = z — z. Parameter 5 describes the velocity 
u of the soliton where u — (1 + k,5 2 )/(1 — k5 2 ). Then, 



A = ± I J x-ut) (72) 



u 



where ± denotes the sign of S. 
A few remarks are in order. 

(i) For iV = —N\ the two expressions of gf 1 ^ given in Eqs.(67) and (68) yield the same 
result, Eq.(70), thus proving the consistency of two expressions in this case. 

(ii) We may obtain an abelian limit by taking a\ = 0,2 = 0, 03 = 1 and g\ = g^ 1 = 
exp(iy/7T(p<j3) . In which case, Eq.(70) reduces to the well-known 1-soliton solution of the 



sine- Gordon equation ||T6 



<p = -A ta n-V A . (73) 

' 7T 



(hi) In the parametrization g 1 1 g2 = exp(i0ajO"i), 1-soliton can be written by 

<h = 2cos~ 1 (-tanh2A) = 2sin _1 — ^— - . (74) 
Y v ; cosh 2 A v ; 
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Note that A changes from ±00 to =Foo as x goes from —00 to 00 so that the soliton number, 
(0(oo) - 0(-oo))/2?r, is ±1. 



M=2; soliton(antisoliton) - soliton(antisoliton) scattering 

For M=2, two possible solutions of Eq.(62) are 

(i) Nl = -N u N\ = -N 2 

(ii) N\ = -N 2 , (75) 

which describe 2-soliton solutions and nonabelian breather solutions respectively. First, we 
consider the case (i). We parametrize Ni,N 2 by 

Ni = % exp(2A! + ri^am , N 2 = i exp(2A 2 + 772)6^ (76) 

where = i5k, A fc = —Kb k z — z/5k and Oj,6j, rji are real constants with normalzation 
didi = bibi = 1. In order to check that the criterion (i) indeed satisfies Eq.(62), we note 
that, for example, Z 2 \Z 22 Z\ 2 = Z\ 2 Z 22 Z 2 \ due to the property that NiNj = exp(4Aj + 2r)i) 
which is proportional to the identity matrix. Therefore, 

(Z 1 )li = (Zli — Z 2l Z 22 ^ Z\ 2 ) 1 = [Zu — Z 2 iZ 22 l Z12) 1 = (Zu — Z 12 Z 22 1 Z 2 i) 1 = {Z 1 ) 11 . 

(77) 

Similar procedure for other components of Z 1 leads to Eq.(62). 

We now calculate g~\ X g 2 for the 2-soliton solution. From Eqs.(66) and (76), we have 

y = (dety) U/i)(i+t>-' id-W)* ) (78) 

where 

detY = -(I + e 4Al+2r?1 )(l + e 4A2+2 " 2 ) - — ^ - (1 - Aq7V 2 )(l - N 2 N ± ) . (79) 
4 (0i + b 2 y 

Thus, gi 1 g 2 can be readily calculated from Eq.(67) or Eq.(68). Either case gives rise to the 
same result in the form; 

g^g 2 = (A + Bia^A + Bicji)- 1 

EE M + MjCTj . (80) 



With the notation, 



A ± = (Ai + ^) ± (A 2 + \ V2 ) , R ee 1-^1 , (81) 
14 



each coefficients are given by 



A = e 2A +[ J R 2 cosh 2 A + -sinh 2 A_ + ^^(aA-l)] 
B\ = iRe 2A+ [Tai sinh(2A 2 + r? 2 ) ± b { sinh(2Ai + 771) + e iik a-b k \ 



(82) 
(83) 



and 

M 



A 2 - B k B k 



A 2 - B k _B k _ 
AR 2 ? 



1 - 



Mi 



n 

A(B% - B 



2 sinh 2 A_ cosh 2 A + + (cosh 2 A + - cosh 2 A_)(l - aA) 
- %e ldk B\B k L 



(84) 



2iR 



A 2 - B k B k 

{Rbi - (P_ + Ra k b k )a t } sinh(2A 2 + r] 2 ) + {Rm + (P_ - Ra k b k )b t } sinh(2A! + 771) 

(85) 



where 



1 t Ft 2 

P± = (R 2 cosh 2 A+ ± sinh 2 A_ t — (aA - 1)) 



(86) 



As in the 1-soliton case, we make the choice of the origin of the coordinate to set parameters 77; 
to zero. Parameters <5j also describe the velocity of solitons. As we show below, if 5 2 = 1/k8\, 
it descibes the soliton - soliton, or antisoliton - antisoliton scattering, whereas if 5 2 = —1/kSi, 
it descibes the soliton - antisoliton scattering in the center of mass frame. In both cases, 
velocities of each solitons are given by u = (1 + K5f)/(1 — K5f) and —u. In the soliton - 
soliton scattering case, R — —1/u , \J — k5\ = tV^ — + u and 



Ai - A 2 



-KX 



=F 



A+ = A! + A 2 = ± 



w 

2J—KUt 



±T . 



(87) 



The upper sign corresponds to the soliton - soliton and the lower sign to the antisoliton 
antisoliton scatterings respectively. In the soliton - antisoliton case, R = —u and 



Ai - A 2 



-KUt 



VT^u 2 



±T 



A+ = A 1 + A 2 = t I = = TX . 



(88) 



where the upper(lower) sign corresponds to the soliton(antisoliton) - antisoliton(soliton) 
scattering. In order to convince the correctness of the solution given by Eqs.(80) - (85), we 
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have checked explicitly that Eqs.(80) - (85), together with 

= 2k{(Y- 1 ) u 5 1 N 1 + 12^2^2 + (Y-'^N, + {Y- l ) 22 5 2 N 2 } , (89) 

indeed satisfy the matrix sine-Gordon equation (12). However, instead of giving cumbersome 
details of the calculation, we present another consistency check. In the abelian limit, where 
we take a\ — b\ — 1, a 2 = a 3 = b 2 = b 3 = and g x = g^ 1 = exp(iy/7iipai), Eqs.(80) -(85) 
gives rise to 

iMi iA{B\-B 1 _) 2iAB\ 



tan(2A/7T9?) 



M 



A 2 - B\B\ " A 2 + (Bl) 2 



and 



2i 1 
B\_ = ±— exp(±2T) sinhXcoshT , A = exp(±2T)(— cosh 2 T - sinh 2 X) 



a 



u 



Using the identity, 



tan 49 



sin4# 



4 sin 9 cos #(cos 2 9 — sin 2 i 



we obtain 



1 



■ tan 



cos 49 (cos 2 9 - sin 2 9) 2 - 4 sin 2 9 cos 2 9 ' 

_ x wsinhX 



tan 49) 



1 



-49) 



tan 



(90) 



(91) 



(92) 



(93) 



2^/ii 2^pn -s/n coshT 

which is precisely the 2-soliton solution of the sine-Gordon theory for the soliton-soliton scat- 



tering with the plus sign and the antisoliton-antisoliton scattering with the minus sign [16 



In order to have a pictorial description of scattering of solitons, we take without loss of 
generality, {aj} = (1,0,0) , {bi} = (cos#, sin 9, 0). Then, 

4 



M 
M l 

M 2 

M 3 
where 



1 - 



u 2 P 2 
2% 



2 sinh 2 X cosh 2 T + (1 - cos #)(cosh 2 T - cosh 2 X) 



uP 2 
2i 



P_{sinh(T + X) - cosflsinh(T - X)} + - (1 - cos 2 9) sinh(T - X) 



-uP_ sin 9 sinh(T — X) + sin 9 sinh(T + X) — cos 9 sin 9 sinh(T — X) 



1 



1 



P± = — cosh 2 T ± sinh 2 X ± -(1 - cos 9)(1 =f -r ) . 

It 2 M 



(94) 



(95) 



The internal motion of solitons may be described most naturally in terms of the parametriza- 
tion g\ X g 2 = exp(— 2iy/7i(piai) where ipi is related to by 



Mo = cos2v / vrM , Mj = — z-j— 7 sin 2y/7r\(p\ 

m 



1,2 



(96) 
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Figures (l)-(3) show M , M\fi , M 2 /i, as an example, for a specific case where 9 = OAn and 
the velocity u = 0.1. From Eq.(7), the potential energy can be written by V = 2Nk,M /-k so 
that Mq depicts the trajectory of soliton - soliton scattering in terms of minus the potential 
energy. M shows that two solitons repulse each other at the origin. It is easy to read the 2tt- 
angle variation A\2^/^T^p\ = 2tt across each bump from these figures which shows clearly that 
they describe the soliton - soliton scattering. Note that the internal direction given by the 
vector {(fij} changes after the collision. Under the spacetime inversion; (T,X) <-> (— T, — X), 
the solution in Eq.(94) possesses symmetry; M <-> M ,Mj <-> — ; i — 1,2. Thus, the 
internal directions of each solitons, specified by the components tpi, become exchanged in 
the process of scattering. This is a characteristic of the scattering of nonabelian solitons. 
The minimum points of M constitute a trajectory of the center of each solitons. At time 
T = — T and T = T , two solitons are located at ±X which satisfies the relation, 



sinhX 2 = V .,2^ V , + 7T2 cosh2T o + £(1 - cos W - ~ 2 )- (97) 



2(cos^ - 1) (cosh 2 T - 1 

2 cosh 2 T - 1 + cos 9 ' a 2 " ' 2" (/ 

Notice that at Tq = 0, two solitons approach closest with 



, 1 1 + cos^ 1 — cos^ 

= ___ + __. (98) 

This shows that the repulsion between two solitons becomes maximum when two vectors 
Oj, bi are aligned in the same direction which is precisely the abelian case where 9 = 0. If 
9 = 7r, Xo takes a minimum value thereby maximizing the nonabelian effect. On the other 
hand, when T becomes large, Eq.(97) can be approximated by 

I 

sinh X 2 w — cosh 2 T . (99) 

The elapsing time for two solitons to bounce back to the separating distance 2X is given by 
2T . Thus, when X becomes large, the elapsing time becomes independent of the angle 9. 



M=2, nonabelian breather 

Now we consider the case (ii); N\ = —N 2 in Eq.(75). We take /ii = — = ^=exp(ia) 
and parameterize N\ and N 2 by 

Ni = iexp(2A)(c k + id k )a k 

N 2 = iexp(2A*)(c k -id k )<r k (100) 

where qq = didt = l 7 and 

2A = 2 V— n[exp(ia)z — exp(— ia)z] = 2\/ — k[(cos a)x + i(sia a)t] 
7 This normalization merely dictates the choice of the origin of coordinates x and t. 
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A straightforward calculation shows that Eq.(62) holds for Ni,N 2 given in Eq.(lOO). We 
could follow a similar procedure as in the case of the two soliton scattering and make use of 
the fact; Z\ 2 = Z 21 , Z\ x = Z 22 and Nf is proportional to the identity matrix. 
From Eqs.(66) and (100), we have, 



Y -i = /, fY) -i / (I + iViiVi t )[l + exp(-2m)]- 1 ^exp(4A*) - \ 

" 1 ' \ -iadi exp(4A) -\ (1 + N\Ni) [1 + exp{2ta)]- 1 



where 



det Y = — [1 + 4(cidi) 2 e 4Kx } + (1 + v 2 )e 2Kx + Ci d t e 2Kx sin(2Kvt). 



(102) 



(103) 



Then, g 1 1 g 2 for the breather solution can be obtained from Eq.(67), or consistently from 
Eq.(68), which we write in the form; 



where 



ih fJ , = (A + B^A + Bla,)- 1 



(104) 



and 
M = 
Mi = 



A = 



2AKx] 



l)e 2Kx - Cl die 2Kx sin2Kvt 



B 



± 



-[1 + 4M0 e 

±iCiv(e Kx sin Kvt + 2c k d k e 3Kx cos Kvt) - 2ive 2Kx e ijk Cjd k . 
±idiv(e Kx cos Kvt + 2c k d k e 3Kx sin Kvt) 



1 - 



2v 2 



(det Y f 



(105) 
(106) 



e 2Kx {l + A(cidij 2 e 4Kx }(l + cA sm2Kvt) + 4c^ i e 4Ka; (c^ i + sin2Kvt) 



\2AKx 



4Kx, 



e Kx sin Kvt + 2c k d k e iKx cos Kvt){2iAdv + 4iv 2 e 2Kx (di - c k d k d)} 



J>Kx 



2 2Kx, 



(det Yf 

+ (e Kx cos Kvt + 2c k d k e 3Kx sin Kvt){2iAdiV + 4iv 2 e 2Kx (c k d k di - a)} 

In addition, a straightforward calculation shows that 
Q^dgi = 2iK,J2( z ~ 1 ) a p 
2i\/—K 



(107) 



detF 



Hi 



ia+2Ai 



(— — + ivc k d k e )(cj + idj)<jj + c.c. 



(108) 



which together with Eqs.(104)-(107) satisfies the matrix sine-Gordon equation (12). 
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For a pictorial description of a nonabelian breather, we choose without loss of generality 
{cj = (1,0,0) , {di} = (cos fl, sin 0,0). Then, 



M n = 1 



2v 



2 r 



e 2Kx ( 1 + 4 cog 2 ^ e 4x^ ( ! + cos ^ sin 2Kwt) + 4 cos fle 4Kx (cos + sin 2tfui ) 



Mi = 



2iu 



(P_e Ka; + 2 cos 2 0P_e 3Ka: - 4v cos sin 2 fle 5 * 1 ) sin Kvt 



Mo = 



+ (cos 6P-e Kx + 2 cos 6P-e 3Kx - 2v sin 2 fle 3 **) cos Kvt 

2iv 
pi 



(sin 26P-e 3Kx + 2v sin #e 3Xx + 2v cos sin 26e 5Kx ) sin lifui 



M 3 = 
where 



+ (P_ sin 0e K * + v sin 2#e 3Xx + 2w sin 29e 5Kx ) cos Kvt 
0, 



P ± = — [1 + 4 cos 2 fle 4X:c ] ± (1 ± t; 2 )e 2X:c ± cos fle 2X:c sin(2i^t). 



(109) 



(110) 



Figures (4)-(6) show M , M 1 /i , M 2 /i for a specific case where = n/6 and t> = 0.1. 
The potential energy profile given in terms of M shows clearly the breathing motion. The 
behavior of ^ in Eq.(96) along the X-direction in the figures of M , M x and M 2 confirms that 
the breather solution is indeed a bound state of soliton and antisoliton. In addition, Mi and 
M 2 shows that the internal direction of the breather also oscillates which is a characteristic 
of a nonabelian breather. Two particular values of 9 are worth to address. If 9 = 0, we 
may follow a similar procedure as in the 2-soliton case, and see that the nonabelian breather 
reduces to the well known sine-Gordon breather, 



y? = ± 



a/? 



tan 



am(Kvt + f ) 
1 cosh(fTx + \n\^2) 



(111) 



For 9 = 7r/2, it is interesting to note that Mo becomes independent of t while Mi and M 2 
are not. This shows that the nonabelian breather at 9 = 7r/2 breathes only internally. That 
is, the internal direction oscillates while the potential energy remains static, i.e. externally 
it becomes completely breatheless. 



4 Backlund Transformation 

The Backlund transformation(BT) is a mapping between two solution surfaces of certain 
differential equations. For example, the sine-Gordon equation dd(f> = sin is invariant under 
the BT 

90 = <90-25sin(^^) 
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B4> = -<90+|sin(^-^) (112) 
o 2 

where 5 is a nonzero real parameter. The integrability of Eq.(112) is the requirement that 4> 
and are both solutions of the sine-Gordon equation. Thus the BT generates a new solution 
from a known one. Moreover, through the Bianchi's permutability theorem, it leads to a 
nonlinear superposition of solutions which gives rise to a new solution by purely algebraic 
means. For example, if a , 0& are two solutions generated by the BT from a known solution 
(f) g with Backhand parameters 8\, 8 2 respectively, then the Bianchi's permutability theorem 
gives a new solution (fi h by 

0h = 5 + 4tan (- —tan ) (113) 

Ox — 6 2 4 

In this section, we show that all these properties generalize to the matrix sine-Gordon 
theory. Recall that the linear equation for the matrix sine-Gordon equation is 

[d + g- l dg - XT]^ g = ; [8+ ~g- x fg]*g = (114) 

where T,T are as in Eq.(20). The BT between two solutions g and / of the matrix sine- 
Gordon equation may be defined in terms of 

*'=A7U< 1 + ^ lf ^ (115) 

where 5 is a real Backhand parameter, tyf and ty g both satisfy the linear equation with 
respect to / and g. On the other hand, if Eq.(115) is combined with the linear equation 
(114) to elliminate then we have an equivalent expression for the BT in terms of / and g, 

f- x df - g-'dg + Sig^Tf , T] = (116) 

and 

gf- 1 ? - Tgf- 1 - 8dgg- l T + 5T Off- 1 = (117) 
Also, the unitarity condition, = 1, requires that 

f- 1 fg-g- 1 ff = 0. (118) 

Eqs. (115) - (118) consitute the BT for the matrix sine-Gordon equation. With g , T , T 
as in Eqs. (20) and (21), Eqs. (116) and (117) in block components are 

fr 1 dfi-gI 1 dg 1 + K5(g 2 1 fi-gi 1 f2) = 

f2~ 1 9f2-92 1 dg2-^(g^f 1 -g^f 2 ) = (119) 
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and 



dhfr 1 -dg 2 92 1 + ](92f2 1 -9ifr 1 ) = 

9 fife 1 - Bgigi 1 - ^fe 1 - gifr 1 ) = 0, (120) 

while Eq.(118) becomes 

f 1 - 1 92 = 9i 1 f2- (121) 

Note that Eqs.(119) and (120) are consistent with the constraint equation (13). We now show 
that the matrix sine-Gordon theory admit also a nonlinear superposition rule of solutions 
which generates a new solution by purely algebraic means. This is given by the permutability 
of the Bianchi diagram. Let a and b are solutions of the matrix sine-Gordon equation 
generated by the BT from a known solution 9 with the Backlund parameters Si and S2 
respectively. Further, let h and h' denote solutions obtained by applications of the BT 
with parameter 82 to a and with parameter 81 to b. Then, the permutability of the Bianchi 
diagram requires h = h . In terms of the BT in Eq.(115), this means that 

- irk atV + r lf 

- XT* AT*< 1 + T r,ft )< 1 + T rlf «»*» < 122 » 

or 

(1 + jh-'fa^l + ja-'fg) = (1 + 6)(1 + ^fg) (123) 

which, when solved for h using the relation TgaT x T~ l = ag^ 1 , gives 

h = giSib' 1 - 8 2 a- 1 )(8ia- 1 - 62b- 1 )- 1 . (124) 

or, in terms of hi and h 2 , 

hi = gi{8ib^ - 8 2 a^){8ia^ - 82K 1 )- 1 

h 2 = g 2 {8ib^ - 8 2 a^){8ia^ - 8 2 b^)- 1 . (125) 

It is easy to check that h is unitary if a, b are unitary. This is the nonlinear superposition rule 
of the matrix sine-Gordon equation which allows one to generate a new solution from a known 
one by purely algebraic means. In the abelian limit, we may take h = expi0^a"3/2 , b = 
expi(ft b (j 3 /2 , a = expi<f> a a 3 /2 , g = expi<f> g a 3 /2, then Eq.(124) reduces precisely to the 
nonlinear superposition rule of the sine-Gordon equation in Eq.(113). 

Finally, we obtain one and two soliton solutions of the theory using the BT. We take the 
trivial solution to be a vacuum given by fi — f 2 — f for a constant SU (2) matrix /. Then, 
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the 1-soliton solution in terms of gi, g 2 is obtained through the BT in Eqs.(119) and (120) 
which, after redefining g u g 2 by g 1 f~ 1 -> g u g 2 f~ 1 -> g 2 , becomes 

g^dgi - re<5(02 1 - g^ 1 ) = 

dgig^ + ^-gi) = 0. (126) 

and the same equation with gi and g 2 interchanged. If we use the parametrization for g\ 
and g 2 , 



f u i\J\ — uu*e l ^\ _( v iy/l — vv*e l6 \ / 107 \ 

9l ~\iy/l^H*e-^ u* ) ' 92 -\i^T=W*e- i6 v* J' [U7) 

Eq.(126) resolves into the component equations; 

u*du - udu* - 2i(l - uu*)d(f) - 2k5(v* - u*) = 

du - K5(e i( * _0 Vl -uu*Vl -vv* - 1 + uv*) = 
m<9m* -M*9w-2i(l-MM*)90 + 2i(M-t;) = 



du - ^(e i(fl_ *Vl - uu*y/l -vv* -1+ uv*) = (128) 

and the same equation with the interchange; u v,(f> 0. In addition, the traceless 
conditon of Eq.(126) requires that u + u* — v + v*. The unitarity condition, Eq.(121), in 
this case requires that u — v* and <fi = 9 + n. Then, equations for u and u* become 

(^-d-Sd)u = 

KO 

(—d + 5d)u = -4 + 2uu* + 2u 2 (129) 

KO 

and the same equation with u and u* interchanged. This equation can be readily integrated 
to give u + u* = 2tanh(2A + rj) ; A = —k5z — z/8 and 

u = tanh(2A + rj) + ic — — \ r, (130) 

v ' cosh(2A + r?) v ' 

where c is an arbitrary constant. In addition, when the solution u is used, Eq.(128) can be 
solved for such that is a constant. In terms of g\ and g 2 , this means that 

9i ~ 92 -y^ 

N = ie 2A+n a k (T k , V=(l + k5 2 )/{\ - k5 2 ) 

A = -k5z-\z = ± /~" (x-Vt) (131) 
5 VT :r V 2 
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where T], dj are arbitrary constants coming from c and with normalization a^ai = 1. This 
agrees precisely with the 1-soliton solution in Sec. 3. 

In order to obtain two soliton solutions, we may apply the nonlinear superposition rule 
to a couple of one soliton solutions obtained by the BT with parameters Si and 82 such that 

ai = a 2 = — , ; 61 = o 2 = — (132) 

2 1 + Nt ' ' 1 2 1 + n 2 K ' 

where Ni, N2 are given in Eq.(131) with repective parameters 8i,t]i and 82^2. Then, from 
Eq.(125) we obtain the 2-soliton solution, 



1-iVi I-JV2/V l-JVa l-JVi 

^ - - 82^) (8^ - 8 2 1 -^Y 1 . (133) 

2 \ 1 1 + N 1 2 l + iV 2 A l + iV 2 2 l + JVi/ 1 ; 

It is now a straightforward but amusing exercise to check that h^ 1 is equal to (A+B^_ak) / det Y 
of Eqs.(79) and (80) while h^ 1 is equal to (A + B k a k )/ det F. 



5 Discussion 



Throughout this Paper, we have analyzed various classical properties of the matrix sine- 
Gordon theory for the coupling constant k < 0. For k > 0, the vacuum structure changes, 
i.e. degenerate vacuua occur at <p — (2?i + l)7r for integer n so that (7i" 1 (72 = — 1, o r for 
gi = —g 2 = f for arbitrary /. This reflects the symmetry of the matrix sine-Gordon theory 
under the exchange; 

9\ «-> 9i , 92 «-> -^2 , « <-> -« • (134) 

Thus, the matrix sine-Gordon theory with k > is identical with the k < case except the 
sign change of This type of symmetry arises due to a specific choice of the potential term 
and can be generalized to other types of nonabelian sine-Gordon theory considered in [|J 
0] 0. In the abelian sine-Gordon case, the potential is given by kcos0 and the symmetry 
becomes the well-known one; 

<-> + 7T , K «-> — K . (135) 

In finding exact solutions, solitons and breathers, we have assumed asymptotic boundary 
conditions <7i,<72 —> 1 as x — > ±00. However, this boundary condition may be relaxed 
to accomodate other types of solutions. For example, if we impose a periodic boundary 
condition, i.e. choose the underlying manifold to be a torus, we can not take a nonlocal 
gauge where A = A = except for the trivial holonomy sector of the flat connection A and 
A. Nontrivial sectors could lead to generalizations of solitons and breathers of the matrix 
sine-Gordon theory characterized by the nontrivial holonomy. Another important classical 
aspect of the matrix sine-Gordon theory which we have not considered is the hamiltonian 
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struture of the theory. As in the case of solitons, we expect a nontrivial generalization of 
the hamiltonian structure of the sine-Gordon theory, especially a nonabelian generalzation 
of the R-matrix. 

The quantum matrix sine-Gordon theory is equally important and provides a nontrivial 
generalization of the quantum abelian sine-Gordon theory. Moreover, it is a natural quan- 
tum field theory framework for the Zamolodchikov's integrable $(2,1) perturbation of Ising 
model. In the parafermion case, the S-matrix obtained by Zamolodchikov using the operator 
algebra has been explained nicely in terms of the complex sine-Gordon theory [|17|. Since the 
matrix sine-Gordon theory carries essentially the same structure as the complex sine-Gordon 
theory, one expect that a similar analysis is possible for this case using a semiclassical WKB 
approximation. Work in this direction is in progress and will appear elsewhere. 
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Figure 1: M in soliton - soliton scattering 
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Figure 2: Mi in soliton - soliton scattering 
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Figure 3: M2 in soliton - soliton scattering 
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Figure 4: Mq of nonabelian breather 
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